1;

function [V,S] = geom_sphere(R)
	S = 4*pi*R^2;
	V = (4*pi/3.0) * R^3;
endfunction

function [V,S] = geom_ell(a,b,c)
	V = (4*pi/3.0) * a *b *c;
	local p=1.6075;
	S = 4*pi*( ( (a*b)^p + (a*c)^p + (b*c) ^p ) / 3.0 )^(1.0/p); 
endfunction
